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(54) Method for synchronizing a receiver 

(57) The invention relates to a method for synchro- 
nizing a receiver (1) with a transmitted code-modulated 
spread spectrum signal. The method uses at least one 
reference code (r(x)) ? which corresponds to a code used 
in the modulation. The frequency shift of the transmitted 
signal and the code phase of the code used in the mod- 
ulation are determined in the method. The method also 
includes at least the following steps: 

a sample vector formation step, in which samples 
(1 01 ) are taken from the received signal to form the 



sample vectors (p k (1), p k (2)...p k (N)) f 

a correlation step {102, 103, 104, 105) for forming 

a correlation function matrix (C x k ) on the basis of 

the sample vectors (p k (1), p k (2)...p k (N)), 

an analysis step (106,107) for forming a coherent 

search matrix (A x k ) on the basis of the transpose 

of the correlation function matrix (C x k ), and 

a determination step (110), in which the frequency 

shift and code phase are determined by using the 

values of elements (a x ^/,y)) of said coherent search 

matrix(A xk ). 
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Description 



[0001] The present invention relates to a method according to the preamble of Claim 1 for synchronizing a receiver 
a receiver according to the preamble of Claim 1 7 and an electronic device according to the preamble of Claim 32 
[0002] The GPS system (Global Positioning System) is a well known positioning system, which comprises over 30 
satellites, of which a maximum of 12 satellites are simultaneously visible to a receiver. These satellites transmit various 
information, such as Ephemeris data and clock data parameters which describe the satellite clock and its relationship 
to GPS time. Normally, the receiver used in positioning determines its location by calculating the transit time of signals 
transmitted simultaneously from several satellites of the positioning system to the receiver. To determine its location 
the receiver must typically receive the signal of at least four visible satellites in order to be able to calculate the location ' 
[0003] Each operational satellite of the GPS system transmits a so-called L1 signal at a carrier frequency of 1575 42 
MHz. This frequency is also denoted 154f 0 , where f 0 =10.23 MHz. In addition, the satellites transmit an L2 signal at a 
earner frequency of 1227.6 MHz, that is 120f 0 . In the satellite, these signals are modulated with at least one pseudo- 
random sequence. Each satellite has a different pseudo-random sequence. As a result of modulation, a code modulated 
wideband signal is formed. This modulation technique makes it possible for the receiver to distinguish between signals 
transmitted by different satellites, despite the fact that the carrier frequencies used in the transmission are essentially 
the same. This modulation technique is called Code Division Multiple Access (CDMA). One of the pseudo-random 
sequences used in each satellite for the modulation of the L1 signal is so-called C/A code (Coarse/Acquisition code) 
which is so-called Gold code. Each GPS satellite transmits a signal using a unique C/A code. The codes are formed 

v1o th v 3 m ? dUl( L 2 ° f tW ° 1 ° 23 " bit binary sec > uences - The first binar y sequence G1 is formed using the polynomial 
X +X +1 , and the second binary sequence G2 is formed by delaying the polynomial X 10 +X 9 +X 8 +X 6 +X 3 +X 2 +1 so that 
each satellite has a different delay. This arrangement enables the generation of different C/A codes with a similar code 
generator. The C/A codes are binary codes, the chipping rate of which in the GPS system is 1 .023 MHz The C/A code 
comprises 1023 chips, which means that the repetition time of the code is 1 ms. The carrier wave of the L1 signal is 
further modulated with navigation information at a bit rate of 50 bit/s. The navigation information comprises information 
about the "health" of the satellite, its orbit and clock data parameters, etc. Each satellite of the GPS system maintains 
a local time by e.g. atomic clocks. 

[0004] During operation, the satellites monitor the condition of their equipment during operation For example sat- 
ellites can use so-called watch-dog functions to detect faults in the equipment and provide notification about them 
Errors and malfunctions can be either momentary or long-term in nature. On the basis of the health information some 
of the errors can be possibly compensated for, or the information transmitted by a malfunctioning satellite can be 
ignored ent.rely. In addition, in a situation in which the signals of more than four satellites can be received information 
received from different satellites can be weighted in different ways on the basis of the health information. Thus it is 
possible to minimize errors that satellites which seem to be unreliable may cause in the measurements 
35 [0005] In order to detect the satellite signals and to identify the satellites, the receiver must perform a synchronization 
operation, in which the receiver searches for the signal of each satellite and tries to synchronize with the signal so that 
the data transmitted with the signal can be received and demodulated. 

[0006] In prior art receivers, the time required for this synchronization operation depends among other things on the 
strength of the received signal. Typically, the weaker the signal to be received, the longer each element of the search 
space (correlation/frequency) must be integrated in orderto detect a possible signal. In prior art GPS receivers which 
are typically designed for outdoor use, synchronization with the satellite signals takes some tens of seconds or a few 
minutes, if the received signal is relatively strong, e.g. in the range of -120 to -130 dBm. However, if positioning is to 
take place indoors or in a place where the received signal is attenuated by buildings or other obstacles in the terrain 
the synchronization time increases substantially. Indoors the signal strength is typically approx -150 dBm whereby 
the time needed for synchronization using a sequential search would be hours. This time is clearly too long for per- 
forming a reasonable positioning. Figure 2a illustrates, by way of example, a CDMA signal transmitted from a satellite 
and correspondingly Fig. 2b illustrates this transmitted signal when it reaches the receiver. The signal to be received 
has been attenuated substantially on the transmission link, and it contains accumulated noise. 
[0007] The positioning receiver must perform synchronization when, for example, the receiver is switched on and 
also in a situation where the receiver has not been able to receive the signal of any satellite for a long time A situation 
like this can easily arise in mobile devices, for example, because the device moves and the antenna of the device is 
not always in an optimal position in relation to the satellites, which weakens the strength of the signal arriving at the 
receiver. In urban areas, buildings also have an effect on the received signal, and in addition, can cause so-called 
multipath propagation in which the transmitted signal arrives at the receiver via several different routes for example 
straight from the satellite (line-of-sight) and also reflected from buildings. As a result of this multipath propagation the 
same signal is received as several signals with different phases. 
[0008] The positioning system has two main functions: 
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1 . calculating the pseudo range of the receiver to different GPS satellites, and 

2. determining the location of the receiver using the calculated pseudo ranges and the satellite position information. 
At any given time the satellite position information can be calculated on the basis of the ephemeris data and time 
correction information received from the satellites. 

5 

[0009] The distances to the satellites are called pseudo ranges, since the time is not known accurately in the receiver 
In that case the determination of location and time is repeated, until a sufficient accuracy of time and a sufficient 
accuracy of location has been achieved. Since the time is not known with absolute accuracy, the location and time 
must be determined by linearizing a set of equations for each new iteration. 
10 [0010] The calculation of the pseudo range can be performed by measuring the average apparent transit times of 
the different satellite signals. After the receiver has synchronized with the received signal, the information transmitted 
in the signal is demodulated. 

[0011] Almost all known GPS receivers use correlation methods for calculating the ranges. The pseudo random 
sequences of different satellites are stored or are generated locally in the positioning receiver. A down conversion is 

is performed on the received signal, whereafter the receiver multiplies the received signal with the stored (or locally 
generated) pseudo random sequence. The signal formed as the result of the multiplication is integrated or low-pass 
filtered, whereby the result indicates whether the received signal contains a signal transmitted by a satellite. The mul- 
tiplication carried out in the receiver is repeated so that each time the phase of the pseudo random sequence saved 
in the receiver is shifted. The correct phase is estimated from the correlation result preferably in such a way that when 

20 the correlation result is the highest, the correct phase has been found. In this way the receiver is correctly synchronized 
with the received signal. 

[0012] After synchronization with the code has been performed, the frequency is fine tuned and phase locking is 
carried out. The correlation result also reveals the information transmitted in the GPS signal, which means that it is a 
demodulated signal. 

25 [001 3] The above mentioned synchronization and frequency adjustment process must be repeated for each satellite 
signal received at the receiver. Therefore, this process consumes a lot of time especially in a situation where the 
received signals are weak. In some prior art receivers, several correlators are used in order to accelerate this process, 
whereby more correlation peaks can be searched for simultaneously. In practical applications, it is not possible to 
accelerate the synchronization and frequency adjustment process very much by simply increasing the number of cor- 

30 relators, because thir number cannot be increased indefinitely. 

[0014] Some prior art GPS receivers use an FFT technique to determine the doppler shift of the received GPS signal 
in connection with ordinary correlators. In these receivers, correlation is used to reduce the bandwidth of the received 
signal to 1 0 kHz-30 kHz. This narrow-band signal is analyzed using FFT algorithms to determine the carrier frequency. 
[0015] A GPS receiver and a method for processing GPS signals are presented in international patent application 

35 WO 97/14057. Principally, the receiver presented in this publication comprises two separate receivers, of which the 
first receiver is intended for use in a situation where the received signal strength is sufficient, and the second receiver 
is intended for use in a situation where the received signal strength is insufficient to perform a sufficiently accurate 
determination of position using the first receiver. In the second receiver, the received signal is digitized and stored in 
a memory means, whereby these stored signals are later processed in a digital signal processing unit. The digital signal 

40 processing unit performs convolution operations on the received, digitized signal. The purpose of these convolution 
operations is to calculate pseudo ranges. The number of code sequences (PM frames) stored in the memory means 
typically ranges from one hundred to thousand, which corresponds to a signal having a length of 100 ms to 1 s. After 
this, a stored code corresponding to the code of a satellite to be examined is retrieved from the memory of the receiver 
to be used in the analysis of the received signal. 

^5 [0016] The Doppler shift is also removed in the receiver. The amount of Doppler shift is determined using either the 
first receiver or on the basis of information received from a base station of the GPS system. After this, consecutive 
frames are coherently summed. A fast Fourier transform is performed on the data set produced as the result of this 
summing. A multiplication is performed on the Fourier transform result, using the complex conjugate of the Fourier 
transform of the reference signal stored in the memory means. An inverse Fourier transformis further performed on 

50 the result of this multiplication, whereby a number of correlation results are received. Thus, in this publication, a cor- 
relation is replaced by a Fourier transform, whereby the number of calculations is reduced. According to the specifi- 
cation, the method accelerates positioning by a factor of 10 to 100 compared with the solutions known at the time of 
filing the application. 

[0017] It is an objective of the present invention to provide a receiver, in which synchronization with a transmitted 
55 signal can be performed substantially faster than in prior art receivers, even if the signals are weak. The invention is 
especially suited for use in positioning receivers, as well as in other receivers, advantageously in CDMA receivers, 
where the receiver must synchroni7e with a spread spectrum signal. The invention is based on the idea that the fre- 
quency analysis properties of convolution and time-to-frequency transformation, such as the fast Fourier transform, 
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are utilized more efficiently, whereby it is possible to examine a greater amount of information simultaneously than in 
pnor art solutions, that employ the Fourier transform. A method according to the present invention is characterized in 
by the which is set forth in the characterizing part of Claim 1 . A receiver according to the present invention is charac- 
terized by that which is set forth in the characterizing part of Claim 17. An electronic device according to the present 
invention is characterized by that which is set forth in the characterizing part of Claim 32. 

[0018] The present invention provides considerable advantages compared with prior art methods and receivers 
Especially m positioning receivers, the method according to the invention enables the receiver to become synchronized 
w.th the received signal faster. In this way positioning can be performed much faster, even in situations where the 
received signal is weak, for example indoors or in urban areas where there are obstacles that substantially attenuate 
the propagation path of the signal. The method according to the invention provides a substantial improvement in the 
synchronization speed when compared with prior art methods. The receiver according to the invention can be imple- 
mented us.ng a relatively small number of components, and the total energy consumption can be maintained at a 
reasonable level, whereby the invention is especially suitable for use in portable devices. Thus the positioning receiver 
can also be implemented in connection with a mobile station. 
is [0019] in the following, the present invention will be described in more detail with reference to the accompanying 
drawings, in which v ' a 

Figure 1 is a simplified block diagram of a receiver, in which the method according to the invention can be applied, 

Figure 2a shows the signal transmitted from the satellite as a simplified diagram, 

Figure 2b shows an example of the simplified signal of Fig. 2a in a receiver, 

Figure 3 shows the correlation step of a method according to a preferred embodiment of the invention as a sim- 

23 plified diagram, 

Figure 4 shows the analysis step of a method according to a preferred embodiment of the invention as a simplified 
diagram, 

30 Figure 5 shows the summing step of a method according to a preferred embodiment of the invention as a simplified 
diagram, 

Figure 6 shows the decision step of a method according to a preferred embodiment of the invention as a simplified 
diagram, ^ 

35 

Figure 7 shows an electronic device according to a preferred embodiment of the invention as a simplified block 
diagram, and 

Figure 8 shows a prior art correlator, which can be used in connection with a method according to a preferred 
40 embodiment of the invention. 



45 



55 



[0020] In the receiver 1 shown in Fig 1 , the received signal is preferably converted to an intermediate frequency in 
he converter block 2. As is well known, at this stage, the signal comprises two components, I and Q, between which 
there is a phase difference of 90». These intermediate frequency converted analogue signal components are digitized 
in digitizing block 3 and conducted to multiplier block 4. In the multiplier block 4, the I and Q components of the digitized 
signal are multiplied with the signal generated by the Numerically Controlled Oscillator (NCO) 5. The numerically con- 
trolled oscillator signal is intended to correct the frequency difference caused by the Doppler shift and the frequency 
difference of the local oscillator (not shown) of the receiver 1 . The numerically controlled oscillator 5 is also used in 
connection w.th the method according to the invention, as will be described later in this description. The signal formed 
by the multiplier block 4 is conducted to the synchronization block 6, where the method according to the invention is 
applied. The synchronization block 6 tries to find the code phase and frequency deviation of the satellite to be used in 
operations performed after the synchronization. This will also be described later in this description. During synchroni- 
zation, the control block 7 controls the scanning block 8 which adjusts the frequency of the numerically controlled 
oscillator 5 as requ.red. The control block 7 causes the first switch 9 to connect either the signal formed by the scanning 
b ock 8 to the numencally controlled oscillator 5 during synchronization or the control signal formed by the code tracking 
block 11 to the numencally controlled oscillator 5 when synchronization has been achieved. The second switch 10 is 
used to control the operation of the code tracking block 11. The code tracking block 11 forms part of a code phase 
locked loop and a carrier phase locked loop (not shown), which are known as such 
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[0021] After switching on the operating voltage, or in a situation where the receiver 1 has been unable to receive the 
GPS satellite signals for a long time, a two-dimensional search step is carried out in the receiver 1 for each satellite 
signal received. The purpose of the two-dimensional search is to determine the carrier frequency and code phase of 
each satellite. The carrier frequency is affected by the Doppler shift resulting from the movement of the satellite and 
5 inaccuracies of the local oscillator of the receiver. The frequency inaccuracy can be rather large, as much as ± 6 kHz, 
in which case the receiver must search a frequency range of approximately 12 kHz with respect to the actual trans- 
mission frequency (L1 = 1575.42 MHz). The receiver 1 does not know the exact code phase, either, whereby the 
receiver must also determine the code phase from 1023 possible code phases. This gives rise to a two-dimensional 
search process, in which frequency deviation in the range of 12 kHz and a code phase out of 1023 different code 

10 phases are searched. In a method according to a preferred embodiment of the invention it is possible to examine a 
frequency range of approximately 500 Hz at a time, whereby the method is repeated 24 times if necessary, to cover 
the whole frequency range of 12 kHz to be examined. Obviously, the values used in this specification only serve as 
examples which clarify the invention, but do not constitute a limitation thereof. The invention can also be applied systems 
other than GPS systems, in which case the frequency values, code phases and the number of codes can vary. 

is [0022] In the following, the operation of a method according to a preferred embodiment of the invention in a receiver 
1 according to Figure 1 will be described. In order to start synchronization, the scanning block 8 sets the frequency of 
the numerically controlled oscillator 5 so that the receiver advantageously receives the lowest frequencies of the fre- 
quency range, which are in this example 1575.414MHz — 1575.4145 MHz. The receiver can also specify the starting 
frequency in such a way that it utilizes previously determined location information and/or almanac information, whereby 

20 the positioning can be further accelerated. Samples of the received signal are saved in sample vector formation block 
12 preferably as complex sample vectors p k (1), p k (2)...p k (N). each of which comprises 1023 samples in this preferred 
embodiment. In this preferred embodiment,, the rate at which samples are stored in the sample vector formation block 
12 is essentially the same as the chipping rate of the chips, which is approx. 1 023 000 samples a second. The sample 
vectors are continuous so that one sample vector continues temporally after another sample vector, that is the time 

25 difference between the last sample of the previous sample vector and the first sample of the next sample vector is 
essentially the same as the time difference between consecutive samples of the sample vector. Thus, the 1 023 samples 
correspond to a signal of 1 ms, which corresponds to a frequency range of 1 kHz of which a part can be utilized. The 
sample vector formation step is denoted with the reference number 101 in the appended Fig. 3. 
[0023] The number of sample vectors is preferably N, where N is advantageously a power of two. In addition, in a 

30 preferred embodiment of the invention the formation of sample vectors is repeated K times, as will be described later 
in this specification. In the following, the subscript k is usedto denote different repetitions. When defining the number 
N of sample vectors p k (1),p k (2)...p k (N). it must be taken into consideration that in the GPS system the signal is mod- 
ulated with information at a bit rate of 50 bits/s as binary phase modulation. r Another factor, which limits the number N 
of sample vectors p k (1), p k (2)...p k (N), is the frequency stability of the receiver's local oscillator. 

35 [0024] In addition to the sample vector formation step, the synchronization method according to the invention also 
comprises a correlation step in which a correlation function matrix is formed. 

[0025] The correlation step can be partly performed already during sampling, or after N sample vectors p k (1 ) ; p k (2)... 
p k (N) have been formed. If the correlation step is performed for example in such a way that a time-to-frequency trans- 
formation, such as a Fast Fourier Transform (FFT), is calculated for each sample vector after it has been saved, the 

40 same time-to-frequency transformer can be used for all N sample vectors p k (1 ), p k (2)...p k (N). If, on the other hand, the 
correlation step is performed after N sample vectors have been saved, a separate time-to-frequency transformer must 
be used for each sample vector, or the time-to-frequency transformations are performed for different sample vectors 
consecutively using the same time-to-frequency transformer. Fig. 3 shows the correlation step of the method, in which 
a correlation function matrix C x k is formed from the sample vectors p k (1 ), p k (2)...p k (N). Although the Fourier transform 

45 is mainly used as an example of a time-to-frequency transformation throughout this description, and the inverse Fourier 
transform is used as an example of an inverse transformation, tha is a frequency-to-time transformation, it is obvious 
that this invention is not limited to these examples only. 

[0026] A discrete Fourier transformation 102, most preferably a fast Fourier transform, FFT is performed on each 
sample vector p k (1), p k (2)...p k (N). 

50 

P k (i) = FFT(p k (i)), where i=1 N (1) 



[0027] This is illustrated by blocks FFT1 , FFT2, .... FFTN in Fig. 3. In practice, the number of values used in calcu- 
lations is preferably 1024, because then the discrete Fourier transformation can be implemented in practical applica- 
tions much more efficiently (with an FFT algorithm) than when 1023 values are used. One way of doing this is to add 
an extra zero as the 1024th element. This has a minor effect on the transformation result. 
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[0028J Preferably, the reference codes r(x) corresponding to the C/A codes of all the satellites of the GPS system 
are stored in the receiver, where x refers to the satellite identifier and is, for example, in the range 1—36. It is not 
necessary to store the reference codes, but they can also be generated in the receiver. In the correlation stage, the 
reference code of the satellite transmitting the signal with which the receiver is being synchronized at a particular time 
is selected or generated. The reference code is temporally reversed. A discrete Fourier transformation 103, most 
preferably a fast Fouriertransform (FFT) is performed on the reversed reference code, which is denoted by f (x) in Fig. 3. 

R(x) = FFT(r{x)) {2) 

[0029] The reversed reference code r (x) and/or its FFT transfom may have been saved in the memory means of the 
receiver in advance, or it is formed from the reference code r(x) in connection with the synchronization process! 
[0030] In the next step of the correlation stage a multiplication 104 is performed between the Fourier transformation 
result P k (i) for each sample vector p k (i) and the Fourier transformation R (x) of the reversed reference code r(x). 



k (i)- R(x).P k (i) (3) 



[0031] An inverse Fourier transformation 105 is performed on the results of these multiplications, whereby the result 
is a cross-correlation of the reference code r(x) and the received signal with all possible integer delays (1023). 

m xk (i) = iFFT(M xk (i)) (4) 

[0032] This result is based on the fact that the Fourier transformation of the convolution of time domain signals 
corresponds to the multiplication of Fourier-transformed signals, that is time domain signals transformed into the fre- 
quency domain. When a reversed reference code is also used, a fast discrete time correlation can be performed using 
a Fouriertransform. Thus, in this preferred example, the cross -correlation result comprises 1 023 elements. The cross- 
correlation results m x k (i) formed from the various sample vectors p k (i) are used to form a correlation function matrix 
C x k , in which the number of rows is the number N of sample vectors. 

[0033] It is obvious that instead of temporal reversal of the reference code, it is possible to form reversed sample 
vectors p k (i) from the sample vectors p k (i), in which case the reference code r(x) is used directly and the reversed 
sample vectors are used in the calculations presented above. In a preferred embodiment, it is not necessary to perform 
either of the reversals mentioned above, but the reference code r(x) and the sample vectors p k (i) can be used as such. 
This is based on using a property of the correlation theorem, which shows that the cross-correlation corr^ , z 2 ) between 
two time discrete functions z A , z 2 can be formed by means of a frequency-to-time transformation of frequency domain 
transformed functions. This can be represented with the following expression: 

40 corr(z 1 , z 2 ) ~ IFFTfZ/ . Z 2 ) (5) 

[0034] Thus, in a preferred embodiment of the invention, it is possible to form the cross-correlation between the 
sample vectors p k (i) and the reference code r(x) by performing a Fourier transform (most suitably a fast Fourier trans- 
form FFT) on the sample vectors p k (i) and the reference code r(x) to form a Fourier transform P k (i) of the sample vector 
p k (i) and a Fourier transform R(x) of the reference code r(x), forming the complex conjugate P* k (i), of the Fourier 
transform P k (i) of the sample vector, multiplying the formed complex conjugate P\(i) with the Fourier transform R(x) 
of the reference code and performing an inverse Fourier transform on the multiplication result. Alternatively, a complex 
conjugate R*(x) can be formed from the Fouriertransform R(x) of the reference code, in which case it is multiplied with 
the Fourier transformation P k (i) of the sample vector, whereafter an inverse Fourier transformation is performed on the 
multiplication result. 

[0035] In this connection, it should be emphasized that the method for calculating the cross-correlation between the 
sample vectors and the reference code r(x) presented in the previous paragraph is a consequence of the basic prop- 
erties of correlation and convolution and the close dependence between them, whereby reversal of a function in the 
^ lime domain is practically equivalent to the formation of a complex conjugate in the frequency domain. This is dealt 
with in more detail e.g. in the publication: "Digital Signal Processing - A Practical Approach", Emmanuel C. Ifeachor 
and Barrie W. Jervis, Addison -Wesley Publishing Company 1 993, ISBN 0-201 -5441 3-X, paragraph 4: "Correlation and 
Convolution", which is referred to in this connection. It should also be mentioned that in general, it is not important for 
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the application of this invention which method is used to obtain the cross-correlation result. 

[0036] The rows of the correlation function matrix C x k formed in the correlation stage represent the cross-correlation 
of the received signal and the reference code with different phase differences at intervals of one millisecond. As a 
formula, the correlation function matrix can be expressed as follows: 
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[0037] In the next phase : that is the analysis phase, the transpose 1 06 of the correlation function matrix C x k is used, 
in which the rows represent the signal samples in the time domain in the same manner as in a prior art correlator. Each 
row corresponds to a certain code phase difference between the received signal and the reference code. A Fourier 
transform 107 is performed on each row of the transpose of the correlation function matrix C x k to form a coherent 
search matrix A x k , whereby a frequency analysis can be performed to determine the actual frequency shift. This is 
illustrated by Fig. 4. 



A x .k=FFT| 



( C I.k) (7) 



25 



[0038] In practical applications, it is not necessary to form a transposed matrix from the correlation function matrix 
in a separate step, but the elements of the stored correlation function matrix C x k are read from memory 1 6 (Fig. 7) in - 
a different direction, preferably by columns. 

30 [0039] The correlation function matrix C x k can also be formed e.g. by using matched filters, in a manner known as 
such. For example, patent specification US-5,373 f 531 discloses a correlator 25 implemented using matched filters 
(Fig. 8), which can be adapted to form the correlation function matrix C x k . The matched filter is formed from three shift 
registers so that the l-component of the received signal is fed to the first shift register 26, the Q-component of the 
received signal is led to the second shift register 27, and the reference code corresponding to the code of the satellite 

35 is fed to the third shift register 28. The length of the shift registers 26, 27, 28 is preferably the same as the length of 
the sample vector, which is 1023 bits in the GPS system. When one sample vector has been transferred to the shift - 
register 26, 27, the correlation is performed as follows. An exclusive NOR operation is performed between the bits of 
the shift register 26 containing the l-component and the bits of the shift register 28 containing the reference code in 
the first reference block 29, and an exclusive NOR operation is performed between the bits of the shift register 27 

40 containing the Q-component and the bits of the shift register 28 containing the reference code in the second reference 
block 30. The result of the operation is bit-specific correlation information: if a bit of the received signal (examined in 
I and Q components) and a corresponding bit of the reference code match, the result is the binary value 1 , and if not, 
the result is the binary value 0. If the signal correlates completely with the reference code, all the bits have the value 
1 . If none of the bits correlate, all the bits have the value 0. The values produced as the result of the correlation are 

^5 summed in the first summing block 31 (l-values) and in the second summing block 32 (Q-values), whereby the result 
is a correlation value m x k (i) having a certain code delay at a certain time. After one comparison, the reference code 
is shifted by one, after which a new correlation result is formed. Each correlation result formed according to the above- 
described method corresponds to one element of a row of the correlation function matrix C x k (1023 in total). The steps 
mentioned above are repeated for each sample vector as required, whereby one row of the correlation function matrix 

so c x k is obtained at each iteration of the calculation. 

[0040] However, as mentioned above, in the GPS system the signal is modulated with a signal of 50 bit/s, which in 
practical applications limits the value of the number N. In this case the number N must be selected preferably so that 
the modulation does not have a substantial effect on the analysis. Furthermore, the optimum value of N is dependent 
on the window function used in the Fourier transformation. For example, if N is chosen to be equal to 32, the the noise 

55 band width is in the order of 30 Hz, which is still a little too large for detecting signals with a strength of the order of 
-150 dBm at the receiver. For this reason, an optional non-coherent summing step, in which the signal-to-noise ratio 
is improved, is performed in the synchronization block 6. 

[0041] In order to implement the non-coherent summing step, the sample vector formation step, the correlation step 
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and the analysis step described above are repeated 108 K times (Fig. 5). The number of repetitions K is preferably 
selected such that the signal-to-noise ratio can be improved sufficiently, in a reasonable time. A coherent search matrix 
^"• k ls formed each time the analysis step is performed, and a non-coherent summing operation is performed on it to 
form a non-coherent search matrix S x . Preferably, the non-coherent search matrix S x is formed as follows Advanta- 
geously, either a magnitude value or some other absolute value, such as the square of the magnitude value, is calculated 
for each of the complex elements a K k (i,j) of each coherent search matrix A, k . The numerical values calculated from 
corresponding elements of each non-coherent search matrix are summed, that is, an addition of the matrices is per- 
formed, which can be expressed with the formula: 
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[0042] In practical applications, a non-coherent search matrix can be formed in at least two ways. In the first alter- 
native implementation, the coherent search matrix formed at each iteration is stored. After the required number of 
repetition, a non-coherent search matrix is formed by summing the equivalent elements according to Formula 8. In this 
alternative implementation, memory is required for storing the elements of all the coherent search matrices. According 
to another alternative implementation, one coherent search matrix is calculated initially, and the values of this matrix 
are copied as the elements of the non-coherent sum matrix. At each iteration a coherent search matrix is formed, and 
the values of that matrix are summed with the corresponding elements of the non-coherent search matrix. In this 
alternative, summing of the equivalent elements is thus performed at each time of iteration. Only one coherent search 
matrix is thus stored, whereby less memory is needed than in the first alternative. 

[0043] After the required number of iterations have been performed, the values of the elements s x (i,j) of the non- 
coherent search matrix S x are examined in the determination step, in an attempt to find 110 a value which exceeds 
the preceding threshold value and is clearly greater than the other values. If such a value like this is found, it represents 
the code phase difference and frequency deviation, because it is probably a signal transmitted by a satellite. If the 
signal is not transmitted by a satellite, but consists of noise or other incidental interference signals, no remarkable 
correlation peaks should arise. The code phase difference is apparent from the row index of the highest value, and 
correspondingly the frequency deviation is represented by the column index. The example of Fig. 6 illustrates a situation 
in which a value, which is clearly larger than the other values, has been found. If, on the other hand, no such value is 
found in the non-coherent search matrix S x , that is, a signal transmitted by the searched for satellite was probably not 
received in the frequency range examined, the frequency range to be examined is changed, and the previously de- 
scribed steps are carried out to form a new non-coherent search matrix. By this method, the whole range of 6 kHz can 
be examined by repeating the steps described above as many times as is required. 

[0044] When required, the previously described steps can be repeated for the whole frequency range to be examined 
and the non-coherent search matrices formed for different iterations or only ranges containing possible peaks can be 
stored before searching for the largest correlation peak. In this way the possibility of false interpretation can be reduced 
for example in such a situation where the threshold value has been set too low and a spurious signal can cause a false 
interpretation. 

[0045] If devices were available, which could identify the code phase and frequency deviation of signals, with a 
strength in the range of -150 dBm, the frequency range would have to be scanned in ranges tens of Hertz in width 
before even an attempt could be made to determ ine the phase and the frequency deviation . Using the present invention ] 
it is possible to scan the frequency range in ranges of hundreds of Hertz in width, that is, much faster than usinq prior 
art receivers. 

[0046] Once the correct frequency deviation and code phase have been determined, the receiver can be set into 
tracking mode. Reception of data is unlikely to be successful for very weak signals and in this case, in a manner known 
as such, it may be necessary to rely on data received via mobile telephone network, for example. In this case distance 
measurement is still possible, but with reduced accuracy. Tracking mode is set by changing the switches (Fig. 1) to 
the second position, whereby received information is fed to code tracking block 11 , in which a feedback connection is 
also provided to enable fine adjustment of the frequency of the numerically controlled oscillator 5. 
[0047] For calculating the location, the receiver performs the signal reception preferably on the basis of signals 
received from at least four satellites. In this case the previously described synchronization process is repeated as 
required on the basis of each satellite signal, whereupon the code of the satellite with which synchronization is per- 
formed is selected as the reference sequence r(x). 
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[0048] Most of the blocks required for the implementation of the method can be implemented for example in a digital 
signal processor (not shown). Either hardware-based solutions or software implementations of the digital signal proc- 
essor can be used to perform FFT transformations. In addition, a control element, preferably a microprocessor or the 
like can be used to control the operation of the receiver. 

5 [0049] An electronic device 24 according to a preferred embodiment of the invention, comprising a mobile station 
and a positioning receiver, is shown in Fig. 7. A first antenna 13 is used for the reception of the signal transmitted by 
the positioning satellites. The received signal is applied to a first radio part 14, in which the signal is down-converted 
and digiii/cd. The first radio part comprises, among other things, the transformation block 2, digitizing block 3 and 
multiplier block 4 of the receiver shown in Fig. 1 . The digitized signal, which at this stage preferably comprises I and 

'0 O components, is applied to the digital signal processing unit 15, in which, among other things, sample vectors are 
formed in the sample vector formation block 1 2. The samples are saved in the first memory means 1 6, which comprises, 
for example, a read/write memory and preferably also read only memory and/or non-volatile read/write memory for 
storing tnc program code of the digital signal processing unit 1 5. In this embodiment, among other things, the functions 
of the synchronization block 6, such as the formation of the correlation function matrix C x k using, for example, Fourier 

'5 transformers FF1 1 . FFT2, .... FFTN and/or matched filters are implemented in signal processing unit 15. Preferably, 
formation ol the coherent search matrix \ k and non-coherent summing step are also performed in the digital signal 
processing unit 15. The digital signal processing unit 15 transmits information about the calculated phase difference 
and the frequency deviation to the processor block 1 7, which comprises, for example, a microprocessor and I/O logic. 
The processor block 17 controls the scanning block 8 and the first switch 9. Preferably, the code tracking block 11 is 

20 at least partly implemented as program instructions of the processor block. The second memory means 1 8 is used as 
the data and program memory of the processor block 1 7. It is obvious that the first memory means 1 6 and the second 
memory means 1 8 can also comprise common memory. Positioning information can be shown to the useron the display 
19. 

[0050] Operations of the mobile station arc also implemented in the application software of the processor block 1 7. 
25 Thus, the display 19 can be used to present telephone call information in a manner known as such. The user can 
control the positioning receiver and the mobile station using the keyboard 20, t. The coding and decoding of audio 
signals is performed by means of a codec 21 . The radio part 22 and the second antenna 23 of the mobile station are 
also shown in Fig. 7 

[0051 ] The present invention is not limited to the previously described embodiments only, but can be modified without 
30 departing from the scope defined by the attached claims. 



Claims 

1. A method for synchronizing a receiver (1) with a transmitted code-modulated spread spectrum signal, in which 
method at least one reference code (r(x)) is used, the reference code corresponding to a code used in the mod- 
ulation ; a frequency shift of the transmitted signal and a code phase of the code used in the modulation is deter- 
mined, characterized in that the method includes the following steps: 

a sample vector formation step, in which samples (101) are taken from a received signal to form the sample 
vectors (p k (1). p k (2)...p k (N)) ; 

a correlation step (102,103,104,105) to form a correlation function matrix (C x k ) on the basis of the sample 
vectors (p k (1). p k (2)...p k {N)), 

an analysis step (106, 107) to form a coherent search matrix (A x k ) on the basis of the transpose of the corre- 
lation function matrix (c x k ). and 

a determination step (110), in which the frequency shift and code phase are determined using values (a xl ^i, 
/)) of elements of said coherent search matrix (A x k ). 

2. A method according to claim 1 . characterized in that in the method a frequency range to be examined is specified, 
the specified frequency range is divided into two or more parts, wherein said sample vector formation step, cor- 
relation step and analysis step are performed on each part, values of elements (a x *(/,/)) of the coherent search 
matrix (A x k ) arc saved, wherein the determination step is performed, and the frequency shift and code phase are 
determined on the basis of a largest value. 

3. A method according to claim 1 or 2. characterized in that a threshold value is specified in the method, wherein 
values of elements (a A k (/ : j)) of the coherent search matrix (A x k ) that exceed said threshold value are used in the 
determination step to determine the frequency shift and the code phase. 
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4. A method according to any of the claims 1 , 2 or 3, characterized in that in the method said sample vector formation 
step, correlation step and analysis step are repeated (108) to form at least two coherent search matrices (A x k ), 
and that a summing step (109) is also performed in the method, in which a non-coherent search matrix (S*)' is 
formed by summing non-coherently the values (a xjk (/ f y)) of the corresponding elements of the coherent search 
matrix (A x k ) formed at each repetition, wherein said'non-coherent search matrix (S x ) is used in said determination 
step to determine the frequency shift and the code phase. 

5. A method according to claim 4, characterized in that in the method a frequency range to be examined is specified, 
the specified frequency range is divided into two or more parts, whereby said sample vector formation step, cor- 
relation step, analysis step and summing step are performed on each part, and values of the elements of the non- 
coherent search matrix (S x ) are stored, wherein the determination step is performed, and the frequency shift and 
code phase are determined on the basts of a largest value. 

6. A method according to claim 4 or 5, characterized in that a threshold value is determined in the method, whereby 
magnitude values of the elements of the non-coherent search matrix (S x ) that exceed said threshold value are 
used in the determination step to determine the frequency shift and the code phase. 

7. A method according to any of the claims 1 to 6, characterized in that the correlation function matrix (C k ) is formed 
using matched filters (26, 27, 28). *' 

8. A method according to any of the claims 1 to 7 T characterized in that 

- in the correlation step, a time-to-frequency transformation (R(x)) of a reversed code corresponding to said 
reference code (r(x)), and a time-to-frequency transformation (Pk(i)) of each sample vector (p k (1), p k (2)...p k 
(N)) is performed, a multiplication between the time-to-frequency transformation (R(x)) of said reversed code 
and the time-to-frequency transformation (Pk(i)) of each sample vector (p k (1), p k (2)...p k (N)) is performed, a 
frequency-to-time transformation is performed on each multiplication result (M x k (i)), and a correlation function 
matrix (C x k ) is formed from the frequency-to-time transformations (m x k (i)) of the multiplication results (M 

(1) ), and that ' xk 

30 ' in the analysis step a coherent search matrix (A x k ) is formed by performing a time-to-frequency transformation 

on the rows of the transpose of said correlation function matrix (C x k ). 

9. A method according to any of the claims 1 to 7, characterized in that 

35 " in tne correlation step, a time-to-frequency transformation (R(x)) of said reference code (r(x)), and a time-to- 

frequency transformation (P k (0) of a reversed vector corresponding to each sample vector (p k (1), p k (2)...p k 
(N)) is performed, a multiplication between the time-to-frequency transformation (R(x)) of said reference code 
(r(x)) and the time-to-frequency transformation (P k (/)) of the reversed vector of each sample vector (p k (1 ), p k 

(2) ...p k (N)) is performed, a frequency-to-time transformation is performed on each multiplication result (M x k 
40 ('))> and a correlation function matrix (c x k ) is formed from the frequency-to-time transformations (m x k (i)) of the 

multiplication results (M x k (i)), and that 

- in the analysis step a coherent search matrix (A x k ) is formed by performing a time-to-frequency transformation 
on the rows of the transpose of said correlation function matrix (C x k ). 

45 10. A method according to any of the claims 1 to 7, characterized in that 

in the correlation step, a time-to-frequency transformation (R(x)) of said reference code (r(x)), and a time-to- 
frequency transformation (P k (/)) of each sample vector (p k (1), pk(2)...pk(N)) is formed, a complex conjugate 
(R*(x)) is formed from the time-to-frequency transformation (R(x)) of the reference code (r(x)), a multiplication 
50 between said complex conjugate and the time-to-frequency transformation (P k (/)) of each sample vector (p k 

Pk( 2 )-Pk( N )) is performed, a frequency-to-time transformation is performed on each multiplication result 
(M x k (i)), and a correlation function matrix (C x k ) is formed from the frequency-to-time transformations (m x k (i)) 
of the multiplication results (M x k (i)), and that' 

- a coherent search matrix (A x k ) is formed in the analysis step by performing a time-to-frequency transformation 
55 on tne r ows of the transpose of said correlation function matrix (c x k ). 

11. A method according to any of the claims 1 to 7, characterized in that 
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in the correlation step, a time-to-frequency transformation (R(x)) of said reference code (r(x)), and a time-to- 
frequency transformation (P k (0) of each sample vector (pk(1), p k (2)...p k (N)) is formed, a complex conjugate 
(P\(i)) 's formed from the time-to-frequency transformation (P k (/))) of each sample vector (p k (1) ; p k (2)...p k 
(N)) ( a multiplication between each said complex conjugate and time-to-frequency transformation (R(x)) of the 
5 reference code (r(x)), is performed, a frequency-to-time transformation is performed on each multiplication 

result (M x k (i)), and a correlation function matrix (C x k ) is formed from the frequency-to-time transformations 
(m K k (i)) of the multiplication results (M x k (i)), and that 

in the analysis step a coherent search matrix (A x k ) is formed by performing a time-to-frequency transformation 
on the rows of the transpose of said correlation function matrix (C x k ). 

w 

12. A method according to any of the claims 8 to 11 , characterized in that said time-to-frequency transformation is a 
Fourier transformation, wherein said frequency-to-time transformation is an inverse Fourier transformation. 



13. A method according to any of the claims 1 to 12, characterized in that the magnitude values of the elements (a x k 
15 (',/)) of the coherent search matrix (A^ k ) are used in said determination step (110). 

14. A method according to claim 13, characterized in that the squared magnitude values of the elements (a x *(/,/)) of 
the coherent search matrix (A x k ) are used in said determination step (110). 

20 15. A method according to any of the claims 1 to 14, in which the location of the receiver (1) is determined and the 
location information is stored, characterized in that in the method a frequency range to be examined is specified, 
and an initial frequency for beginning the examination of the frequency shift is selected from said frequency range 
by using the location information stored in the receiver (1). 

25 16. A method according to any of the claims 1 to 1 5, in which method a spread spectrum modulated signals transmitted 
by two or more satellites are received in the receiver (1), characterized in that in the method said sample vector, 
formation step, correlation step, analysis step and determination step is performed on the signal transmitted by 
each of said satellites. 



30 17. A receiver (1), which comprises at least a synchronization means for synchronizing the receiver (1) with a trans- 
mitted code-modulated spread spectrum signal, the receiver (1 ) having means (1 6) for using at least one reference 
code (r(x)) in connection with synchronization, which reference code (r(x)) corresponds to a code used in the 
modulation, and means (15) for determining a frequency shift and a code phase of the code used in the modulation, 
characterized in that the receiver (1) also comprises; 

35 

sample vector formation means (12, 16) for forming sample vectors (p k (1), p k (2)...p k (N)) from the received 
signal, 

correlation means (6, FFT1 , FFT2...FFTN) for forming a correlation function matrix (C x k ) on the basis of the 
sample vectors (p k (1), p k (2)...p k (N)), 
40 - analysis means (6, 12) for forming a coherent search matrix (A x k ) on the basis of the correlation function 

matrix (C x k ), and 

determination means (7, 1 7) for determining the frequency shift and code phase using the values of elements 
( a x. *(''/)) °f said coherent search matrix (A x k ). 



45 18. A receiver (1) according to claim 17, characterized in that it comprises means (17) for specifying the frequency 
range to be examined, means (5, 8, 9) for dividing the specified frequency range into two or more parts, wherein 
the formation of sample vectors, the formation of a correlation function matrix (C x k ) and Ihe formation of a coherent 
search matrix (A x k ) are arranged to be performed on each part, means (1 6) for saving the values of the elements 
(a x k (ij)) of the coherent search matrix (A x k ), and means (15) for determining the frequency shift and code phase 

so on the basis of a largest value. 



19. A receiver (1) according to claim 17 or 18, characterized in that it comprises means (15) for specifying a threshold 
value, and means (1 5) for comparing said threshold value and the values of said coherent search matrix (A x k ) for 
determining the frequency shift and code phase; 

55 

20. A receiver (1) according to any of the claims 17, 18 or 19, characterized in that the formation of sample vectors, 
the formation of a correlation function matrix (G xk ) and the formation of a coherent search matrix (A x k ) are arranged 
to be repeated for forming a at least two coherent search matrices (A x k ), and that the receiver (1) also comprises 
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summing means for forming a non-coherent search matrix (S x ) by summing non-coherently the values (a xk (i,j)) of 
the corresponding elements of the coherent search matrix (A x k ) formed at each repetition, wherein the*receiver 
(1 ) also comprises means for determining the frequency shift and code phase by using said non-coherent search 
matrix [S x ). 

21. A receiver (1 ) according to claim 1 7, characterized in that it comprises means (1 7) for specifying the frequency 
range to be examined, means (5, 8, 9) for dividing the specified frequency range into two or more parts, wherein 
tnc formation of said sample vectors, the formation of a correlation function matrix (C x k ), theformation of a coherent 
search matrix (A x k ) and the summing of the coherent search matrix (A x k ) to the non-coherent search matrix (S x ) 
arc arranged to be performed on each part, means (16) for saving the values of the elements of the non-coherent 
search matrix (S x ), and means (15) for determining the frequency shift and code phase on the basis of the largest 
value 

22. A rcccivet (1) according to claim 21, characterized in that it comprises means (15) for specifying the threshold 
value, ano means (1 5) for comparing said threshold value and the values of said non-coherent search matrix (S x ) 
tor determining the frequency shift and code phase. 

23. A receiver (1) according to any of the claims 17 to 22, characterized in that it comprises at least one matched 
ftller (26. 27. 28) for forming a correlation function matrix (C x k ). 

24. A receiver (1) according to any of the claims 17 to 23, characterized in that 

the correlation means comprise means (1 6) for forming a time-to-frequency transformation (R(x)) of a reversed 
code corresponding to said reference code (r(x)), means for forming a time-to-frequency transformation (P k 
(i)) of each sample vector (p k (_1), p k (2)...p k (N)), means (15) for performing a multiplication between the time- 
to-frequency transformation (R(x)) of said reversed code and the time-to-frequency transformation (Pk(i)) of 
each sample vector (p k (1), p k (2)...p k (N)) ; means (15) for performing a frequency-to-time transformation on 
each multiplication result (M xk (i)), and means (15) for forming a correlation function matrix (C x k ) from the 
frequency-to-time transformations (m x k (i)) of the multiplication results (M x k (i)), and that 
the analysis means comprise means (15) for forming a coherent search matrix(A x k ) by performing a time-to- 
frequency transformation on the rows of the transpose of said correlation function matrix (C x k ). 
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25. A receiver (1) according to any of the claims 17 to 23, characterized in that 

35 " the correlation means comprise means (16) for forming a time-to-frequency transformation (R(x)) of said ref- 

erence code (r(x)). means for forming a time-to-frequency transformation (P k (/)) of a reversed vector corre- 
sponding to each sample vector (p k (l), p k (2)...p k (N)) : means (15) for performing a multiplication between the 
time-to-frequency transformation (R(x)) of said reference code (r(x)) and the time-to-frequency transformation 
(P k (/)) of the reversed vector of each sample vector (p k (1 ), p k (2)...p k (N)), means (1 5) for performing a frequen- 
cy-to-time transformation on each multiplication result (M x k (i)), and means (15) for forming a correlation func- 
tion matrix (C x k ) from the frequency-to-time transformations (m x k (i)) of the multiplication results (IVL J\)) and 
that X,K 

- the analysis means comprise means (15) for forming a coherent search matrixQ by performing a time-to- 
frequency transformation on the rows of the transpose of said correlation function matrix (). 

26. A receiver (1) according to any of the claims 17 to 23, characterized in that 

the correlation means comprise means (1 6) for forming a complex conjugate (R*(x)) of the time-to-frequency 
transformation (R(x)) of the reference code (r(x)), means for forming a time-to-frequency transformation () of 
each sample vector (p k (1) ; p k (2)...p k (N)) ? means (15) for performing a multiplication between the complex 
conjugate (R*(x)) of the time-to frequency transformation of the reference code (r(x)) and the time-to-frequency 
transformation () of each sample vector (p k (1) ; p k (2)...p k (N)), means (15) for performing a frequency-to-time 
transformation on each multiplication result (M x k (i)) : and means (15) for forming a correlation function matrix 
() from the frequency-to-time transformations (m x k (i)) of the multiplication results (M x k (i)), and that 
the analysis means comprise means (15) for forming a coherent search matrix () by performing a time-to- 
frequency transformation on the rows of the transpose of said correlation function matrix (). 

27. A receiver (1) according to any of the claims 17 . to 23, characterized in that 
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the correlation means comprise means (16) for forming a complex conjugate () of the time-to-frequency trans- 
formation ()of each sample vector (p K (1),p k (2)„.p k (N)), means for forming a time-to-frequency transformation 
(R(x)) of the reference code (n») : means (15) for performing a multiplication between the time-to-frequency 
transformation (R(x)) of the reference code (r(x)) and the complex conjugate () of the time-to-frequency trans- 
formation of each sample vector (p k (1) : p k (2)...p k (N)), means (15) for performing a frequency-to-time trans- 
formation on each multiplication result (M x k (i)), and means (15) for forming a correlation function matrix () 
from the frequency-to-time transformations (m x k (i)) of the multiplication results (M x k (i)), and that 
the analysis means comprise means (15) for forming a coherent search matrix () by performing a time-to- 
frequency transformation on the rows of the transpose of said correlation function matrix (). 

28. A receiver (1 ) according to any of the claims 24 to 27, characterized in that said time-to-frequency transformation 
is a Fourier transformation, wherein said frequency-to-time transformation is an inverse Fourier transformation. 

29. A receiver (1) according to any of the claims 1 7 to 28, characterized in that magnitude values of the elements () 
'5 of the coherent search matrix () are used in the determination means. 

30. A receiver (1) according to any of the claims 17 to 29, characterized in that squared magnitude values of the 
elements () of the coherent search matrix ()are used in the determination means. 

so 31. A receiver (1) according to any of the claims 17 to 30, which comprises means (14, 15, 17) for determining the 
location of the receiver (1) and means (18) for storing the location information, characterized in that the receiver 
(1) also comprises means (17) for specifying a frequency range to be examined, and means (17, 22, 23) for se- 
lecting, on the basis of the location information stored in the receiver (1), an initial frequency of said frequency 
range for beginning the examination of the frequency shift. 

25 

32. An electronic device (24), which comprises at least a positioning receiver (1 ), synchronization means for synchro- 
nizing the positioning receiver (1) to a transmitted code-modulated spread spectrum signal, the positioning receiver 
(1) having means (16) for using at least one reference code (r(x)) in connection with the synchronization, the 
reference code (r(x)) corresponding to a code used in the modulation , and the electronic device comprising means 

30 (15) for determining the frequency shift of the transmitted signal and the code phase of the code used in the . 

modulation, characterized in that the electronic device (24) also comprises: 

sample vector formation means (12, 16) for forming sample vectors (p k (1), p k (2)...p k (N)) from the received 
signal, 

35 - correlation means (6. FFT1, FFT2...FFTN) for forming a correlation function matrix(C x k ) on the basis of the 

sample vectors (p k (1), p k (2)...p k (N)), 

analysis means (6, 12) for forming a coherent search matrix (A x k ) on the basis of the correlation function 
matrix (C x k ), and 

determination means (7, 17) for determining the frequency shift and code phase by using values of the ele- 
40 ments (a xk (i',j)) of the coherent search matrix (A„ k ). 

33. An electronic device (24) according to claim 32, comprising means (14, 15, 17) for determining the location of the 
electronic device (24) and means (18) for storing the location information, characterized in that the electronic 
device (24) also comprises means (17) for specifying the frequency range to be examined, and means (17, 22, 

45 23) for selecting, on the basis of the location information stored in the receiver (1), an initial frequency of said 

frequency range for beginning the examination of the frequency shift. 

34. An electronic device (24) according to claim 32 or 33, characterized in that it also comprises means for performing 
mobile telephone operations. 

50 

35. An electronic device (24) according to claim 34, characterized in that it comprises means (22, 23) for establishing 
a data transfer connection to a mobile communication network, wherein the means for determining the location of 
the electronic device (24) comprise means (22 : 23) for retrieving information for use in the location determination 
from the mobile communication network. 
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values of elements {a x f ^i,j)) of said coherent search 

matrix(A xk ). 
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